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There are at least two methods to calculate Bk with staggered fermions: one is the two spin trace formalism 
and the other is the one spin trace formalism. We have performed numerical simulations on a 16 3 x 40 lattice m 
full QCD with P = 5.7 and a dynamical quark mass 0.01 in lattice units. We try various sources to select only 
the pseudo-Goldstone bosons and compare the various results. 



1. INTRODUCTION 

The lattice simulation of weak matrix elements 
has been one of the most significant contribution 
to Standard Model phenomenology. Especially, 
the knowledge of Bk which describes K° — K° 
mixing through the AS = 2 electro-weak Hamil- 
tonian is important in order to determine pre- 
cisely the parameters of the Standard Model from 
experimental data. 

There are two methods to transcribe a contin- 
uum weak matrix element (eg Bk) to the lattice 
with staggered fermions [1,2]: one is the one spin 
trace formalism and the other is the two spin trace 
formalism. By writing the four-fermion operators 
as a product of operators bilinear in the fermion 
fields, in the one spin trace formalism each exter- 
nal hadron is contracted with both bilinears of the 
four-fermion operators simultaneously, whereas 
in the two spin trace formalism each external 
hadron is contracted with only one of the bilinears 
[1,2]. The problem is that the operators in both 
formalisms are different on the lattice whereas the 
operators in both formalisms are identical in the 
continuum by Fierz transformation. There have 
been a number of attempts to evaluate Bk on 
the lattice with staggered fermions using both a 
Landau gauge operator and a gauge-invariant op- 
erator in the two spin trace formalism [3]. 

Here we will explain the one spin trace for- 
malism and compare the numerical results be- 
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tween both formalisms. Furthermore we try an 
improved numerical technique (cubic wall source) 
in order to select only the pseudo-Goldstone bo- 
son signal exclusively. The results of the cubic 
wall source are compared with the conventional 
method (Kilcup and Sharpe wall source [4]). 

2. ONE SPIN TRACE FORMALISM 

We follow the notation for the four-fermion op- 
erators of Ref. [1,5]. In the continuum, Bk is 
defined as 

B _ (K° \s 7 ,(l- 75 )ds 7 ,(l- 7s )d\ K°) 
K l(K° \s 7 , 75 d \0)(0 \s 7 , 75 d\ K°) . {) 

In the two spin trace formalism, the four-fermion 
operator in the numerator in Eq. (1) is tran- 
scribed to the lattice [1,5] as a sum of four terms: 

= (V x P)H» a + (V x Pf™ 

+ (A x P)H« + (A x P)Z« (2) 

where V (or A) represents the vector (or axial) 
spin structure, P represents the pseudoscalar-like 
flavor structure and the subscript ab; ba (or aa; bb) 
represents the color indices of the quark fields [1]. 
The expression in Eq. (2) has the same chiral 
behavior in the limit of vanishing quark mass as 
does the continuum operator (this is true even 
for the analytic parts of each term in Eq. (2) sep- 
arately) [5,1]. In addition, O^tjI 

preserves the 

same leading logarithmic behavior as the contin- 
uum AS = 2 operator [6,1]. 
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Figure 1. Comparison of naive Bk in both one 
spin trace and two spin trace formalisms calcu- 
lated with Kilcup and Sharpe wall source. 
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Figure 2. Comparison of tadpole-improved renor- 
malized Bk at Tr/a scale in both one spin trace 
and two spin trace formalisms calculated with 
Kilcup and Sharpe wall source. 



In the one spin trace formalism, the four- 
fermion operator of the numerator in Eq. (1) is 
transcribed to the lattice as follows: 



fr\La,tt 
U 1TR 



{V x P)E".f a + (V x P)£* 
+{A x P)%« + (A x P)£* 
_i_ r> 1T . R (V\ 

~ chiral partner * \ / 

In contrast with the two spin trace formalism, the 
individual terms in Eq. (3) do not possess the 
same chiral behavior as the continuum AS = 2 
operator [1]. We must add 0^ R al partner in order 
to preserve the correct continuum chiral behavior 
[1]. By imposing the correct chiral behavior on 
®itr> we determine the chiral partner operator 
[1] as follows: 



n iTR 

chiral partner 



= (^S)™+(^x5)™ 
-{A x S)Z& + U x S)££ . (4) 



This forces the resulting operator to respect the 
continuum chiral behavior. The next question is 
whether the additional chiral partner operators 
still have the continuum leading logarithmic be- 
havior. By choosing basis operators which belong 
to the identity representation with respect to the 
90° axial rotation group (a subgroup of the ex- 
act ?7a(1) symmetry group), we find an eigen- 
operator (Eq. (3) and Eq. (4)) which possesses 
the same chiral behavior and the same leading 
logarithmic behavior as the continuum AS = 2 
operator [1]. 



The Fierz transform of Ofy% in Eq. (3) is dif- 
ferent from OfyQ in Eq. (2) by the following 
terms [1]: 



Oitr - Oar* = {{V + A)x(S- T)) 2TR . (5) 

The contraction of ((V+A)x(S-T)) 2TR operator 
with external pseudo-Goldstone Kaons vanishes 
as a — > 0, since the pseudo-Goldstone Kaon has 
the pseudoscalar-like flavor matrix and the flavor 
trace in the contraction is zero. The above argu- 
ment of vanishing flavor trace is not true for finite 
lattice spacing a ^ since flavor symmetry vio- 
lation (flavor mixing) is present for a ^ 0. This 
means not only that the difference of the numer- 
ical results in both formalisms shows how close 
we are to the continuum but also that as a — > 
the continuum Fierz transformation property of 
AS = 2 operator is recovered [1]. 

The numerical simulation of Bk in both for- 
malisms has been calculated using the Columbia 
16 Giga-flop parallel processor. Configurations 
of 16 3 x 40 size are updated using R-algorithm 
with (3 = 5.7 (1/a = 2 GeV) for two dynami- 
cal flavors with m q a = 0.01. For the Bk mea- 
surement, the configuration is doubled (16 3 x 80) 
along time direction and the distance between 
the two wall sources is 36. The quark mass 
pairs for the Kaon are (0.01,0.01), (0.02,0.02), 
(0.03,0.03), (0.004,0.01), (0.004,0.02), (0.01,0.03) 
and (0.004,0.05). Results calculated with a Kil- 
cup and Sharpe source [4] in both formalisms 
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Figure 3. Naive Bk with wrong (scalar-like) fla- 
vor structure and valence quark mass 0.02 calcu- 
lated with Kilcup and Sharpe wall source in two 
spin trace formalism, which are supposed to van- 
ish as a — > 0. 
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Figure 4. Naive Bk with wrong (scalar-like) fla- 
vor structure and valence quark mass 0.02 calcu- 
lated with cubic wall source in two spin trace for- 
malism, which are supposed to vanish as a — > 0. 



are compared in Figures 1 and 2. The tadpole- 
improved renormalization for Bk are explained 
in detail in Ref. [1]. From Figures 1 and 2, we 
conclude that the tadpole-improved renormalized 
Bk data in both formalisms are in good agree- 
ment with each other. 

3. CUBIC SOURCE METHOD 

For hadron spectrum measurements, the sink 
possesses the same symmetry as a specific 
hadronic state. In contrast to the hadron spec- 
trum measurements, the weak matrix element 
measurements require the symmetry properties of 
the wall source to determine the specific hadronic 
state since the electro- weak effective Hamiltonian 
does not select any particular hadronic state by 
itself. There is a technique (Kilcup and Sharpe 
wall source) which selects only pseudo-Goldstone 
states out of many pion states [4], but may pos- 
sibly also produce p meson states. 

For the cubic source method [7], we implement 
the GRF (Geometric Rest Frame) group on the 
wall source in such a way that for each quark 
flavor, 8 propagators are obtained for sources 
that are non-zero for each corner of the spa- 
tial cubes making up the wall. Specific com- 
binations of these 8 propagators are chosen to 
project on the pseudo-Goldstone mode [7]. We 
have tried these two kinds of wall sources in 



our numerical simulation of Bk- The manifest 
improvement given by the cubic source can be 
seen by examining the wrong flavor channel of 
Br {{V + A) x S) which is supposed to vanish 
in the continuum limit a = 0, as in Figures 3 
and 4. Our results, obtained on 63 configura- 
tions, are naive B 2 K TR = 0.66(6), naive B K TR = 

0. 76.11), tadpole-improved -B^ Tfl (f ) = 0.66(5), 
and tadpole-improved -B^ Tfl (f ) = 0.67(10). 
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